Thermal noise influences fluid flow in thin films during spinodal dewetting 
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Experiments on dewetting thin polymer films confirm the theoretical prediction that thermal noise 
can strongly influence characteristic time-scales of fluid flow and cause coarsening of typical length 
scales. Comparing the experiments with deterministic simulations, we show that the Navier- Stokes 
equation has to be extended by a conserved bulk noise term to accomplish the observed spectrum 
of capillary waves. Due to thermal fluctuations the spectrum changes from an exponential to a 
power law decay for large wavevectors. Also the time evolution of the typical wavevector of unstable 
perturbations exhibits noise induced coarsening that is absent in deterministic hydrodynamic flow. 



PACS numbers: 47.61.-k, 47.15.gm, 68.08.Bc, 68.15.+e 

With the advent of nanofluidics in the last years it 
became evident that thermal noise may play an impor- 
tant role in all hydrodynamic processes occuring at free 
interfaces on small scales. Although in bulk fluids hy- 
drodynamic Navier- Stokes equations are proven to be 
valid down to the nanometer scale, in free interface flow 
stochastic forces induced by molecular motion can sig- 
nificantly alter the behavior even on a micrometer scale. 
Moseler and Landman, for instance, found that the deter- 
ministic lubrication approximation for axial-symmetric 
free boundary flow is not applicable for the description 
of nanoscopic cylindrical jets Q]. The lack of thermally 
triggered fluctuations in the classical hydrodynamic con- 
tinuum modelling was identified as the most likely source 
for deviations of Navier- Stokes equation from molecular 
dynamics simulations. They derived a stochastic differ- 
ential equation that includes thermal noise, whose influ- 
ence on the dynamics increases as the radius of the nano- 
jet becomes smaller, leading finally to the emergence of 
symmetric double cone neck shapes during the breakup 
instead of a long thread solution as expected in the ab- 
sence of noise. In Ref. [2] path integral methods were ap- 
plied to confirm that thermal noise induces indeed qual- 
itative changes in the breakup of a liquid nanometer jet: 
Thermal fluctuations speed up the dynamics and make 
surface tension an irrelevant force for the breakup. Very 
recently, the importance of thermal noise for drop for- 
mation was observed in a colloidal dispersion with an 
ultra- low surface tension [3]. 

Dewetting of thin films is another technologically im- 
portant free surface flow problem where a lubrication ap- 
proximation is used widely. So far, thin film flow has been 
studied by hydrodynamic Navier- Stokes equations with- 
out considering thermal noise, e.g., [4]. With decreas- 
ing film thickness, however, it is expected that thermal 
noise gains significance In Refs. [1, 0, thermal 
fluctuations have been taken into account and a stochas- 
tic version of the thin film equation was derived based 



on the lubrication approximation for stochastic hydro- 
dynamic equations 0, 0]. The magnitude of the added 
stochastic forces increases as the thickness of the film de- 
creases. Recent analytical and numerical studies of the 
stochastic thin film evolution indicate that thermal noise 
increases characteristic timescales of the dewetting pro- 
cess [6] and changes qualitatively the time evolution of 
the film thickness h{r^t) (with r = (x,?/)), i.e., its power 
spectrum C{q^t) [5]. Here, we report that these predic- 
tions of noisy hydrodynamics can be confirmed experi- 
mentally by in situ atomic force microscopy (AFM) of 
dewetting thin polymer films. 

For the experiments, we used polystyrene (PS) of 
2.05 kg/mol molecular weight, prepared on a silicon wafer 
with a 191 nm thick silicon oxide layer. By spin casting a 
toluene solution of PS we prepared films with thicknesses 
of (3.9 ± 0.5) nm and (4.6 ± 0.5) nm (as measured by el- 
lipsometry at different sites). These films are unstable 
and dewet spinodally. Above the glass transition tem- 
perature, the PS film is liquid and capillary waves can 
be observed Below a characteristic wave number 

their amplitudes grow exponentially in time, eventually 
leading to holes in the film [12, 13]. Dewetting then pro- 
ceeds by the growth of holes and their coalescence. By 
tapping mode^^ AFM, the entire dewetting process can 
be monitored in situ, c.f. Fig. [TJ 

In addition to the experiments, we use determinis- 
tic simulation data of the spinodal dewetting process 
from Ref. fl3]. In these simulations, the experimen- 
tal system parameters film thickness, viscosity, and ef- 
fective interface potential ^{h) were used. We gain 
quantitative information about the influence of thermal 
noise to the dewetting process by analyzing the variance 
<j^(t) = — h of the film height h{r^t) and the vari- 
ance k'^{t) (V/i)V(27rcr2(t)) of the local slope of/i(r,t). 
The quantity k'^ {t) gives information about the preferred 
wavevector within the film surface in the early stage of 
dewetting before holes are formed, cf. Fig. [H In this 
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FIG. 1: Dewetting of a 3.9 nm polystyrene film as monitored 
by in situ AFM at 53°C. 



early stage regime, the overbars represent an integration 
over all positions of the image. During the later stage, 
however, when holes appear and grow, the linearization 
of the thin film equation (TFEq) [5], which we use to 
analyse the data, is in general not valid, but may still be 
used in between the holes, where the film heights are close 
to the prepared film thickness h^. We therefore analyse 
in the later stage of dewetting the film only at regions 
where /i(r, t) ^ ho and ignore the other parts. Note that 
this spatially selective data analysis is not possible by 
using correlation function or Fourier space techniques. 

A main problem in the data analysis is the finite size 
of the observation window and the emergence of grow- 
ing holes in the thin film which limits the film regions 
where linear dynamics can be studied. The variances 
cr^ and k'^ (which are moments of the spectrum) cannot 
be determined accurately from the Fourier transform of 
the AFM images. The measured spectrum is convoluted 
with the Fourier transform of the observation window, 
which decays as ~ q~'^ for a sharp cutoff. Thus, the 
short- wavelength behavior of the spectrum is in our case 
dominated by the edges of the AFM image and not by 
the thermal fluctuations. Moreover, the exponentially 
damped spectrum C ~ e"^^!^*^^^!, characteristic for the 
deterministic dynamics, is masked by this convolution 
yielding a spectrum similar to the algebraic noise-induced 
Ccw{q) ^ [5]. Since the 'bare' spectrum C{q) is 
experimentally not accessible, one has to determine vari- 
ances cr^ and k'^ from the real space images. Whereas 
determining is straightforward, calculating k'^ from 
pixelized images using numerical differentiation is too in- 
accurate. 

Here, we apply an a naly sis technique based on 
Minkowski functionals Mj^ [ij] . We threshold the images 
at a certain film height h and measure the area Mo(/i), 
boundary length Mi(/i), and Euler number M2{h) as 
function of h. In the linear regime of spinodal dewet- 
ting, the fluctuations in these quantities should follow 
a Gaussian distribution which was previously tested in 
Ref. [13]. Then, the mean values of Mi,{h) depend only 
on the variances and /c^, which can be determined 
as fit parameters to the functions Mjy{h). In Fig. [2] we 
show the resulting temporal evolution of cr^(t) and of 
k'^{t) in the linear regime. Besides providing accurate 
estimations of these values, this morphometric analysis 



allows to test whether the analysed film region is still 
in the linear regime. As soon as the non-linear dynam- 
ics sets in, the Gaussian model fails to fit the measured 
functions Mj^{h). We compare the data from two experi- 
ments at two spatial resolution lengths (exp. 1: pixel size 
rp ^ 6 nm temperature T ^ 53^C; exp. 2: ^ 40 nm, 
T ^ 65^C) with data from deterministic simulations 
at similar physical conditions and find substantial differ- 
ences. In the following, we explain these differences by 
introducing thermal noise to the description of thin film 
flow. 

In both experiments the short chain length PS film 
can be described as an incompressible Newtonian liquid 
(with a shear viscosity on the order of 77 ~ 10^ Ns/m^) 
on an infinite, flat solid substrate. Thus, the flow can be 
modelled by a stochastic Navier-Stokes equation with an 
additional random stress fluctuation tensor representing 
the effect of thermal molecular motion 0, 0] . We assume 
no-slip boundary conditions between liquid and substrate 
and full-slip at the liquid/air interface z = /i(r, t). At the 
latter, the normal stress is balanced by the surface ten- 
sion 7 and the disjoining pressure — <^>^ For a smooth thin 
film, where the characteristic film height ho ^ 4 nm is 
much smaller than the length scale 100 nm over which 
the film thickness varies laterally, a long wavelength ex- 
pansion yields the following stochastic TFEq 0, IH, 0] 

(1) 

with a single multiplicative conserved noise vector N(r, t) 
obeying (N(r,t)) = and (7V,(r, t) 7V^•(r^ t^) = S^j S{t - 
r') S{t — t'), where (•) denotes the ensemble average over 
realizations of the noise. The characteristic lateral length 
scale for a thin film of /iq = 4 nm is given by the spinodal 
wavelength 27T/qo = ^/-S7T'^J/^"{ho) = ^h\^Ji^^^ijA, 
where the Hamaker constant A determines the effec- 
tive interface potential ^(h) = — i2^h'^ • have 
^ ^ 400 nm with the values A 2 ^ 10"^° Nm and 
7 ^ 3-10~^ N/m determined experimentally for the same 
system [10]. The experimental parameters lead to dimen- 
sionless amplitudes of the noise ^^Jl^^T^ ^ 4 • 10~^ and 

2 • 10~^ for exps. 1 and 2, respectively. In the following 
we show that even this small current induced by ther- 
mal fluctuations can substantially influence the dewet- 
ting process. 

Comparing the experiments with the numerical solu- 
tion of Eq. ([1]) with T = 0, we indeed find excellent 
agreement of the spatial structures Nevertheless, 
the time scales do not match. This can be illustrated by 
comparing the time dependence of cr^(t) and k'^{t) for the 
AFM data and the deterministic thin film simulations, 
see Fig. O The most significant deviation is observed 
for k'^{t) which varies in time for the experiments, but 
stays constant for the deterministic simulations. This 
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FIG. 2: Analysis of AFM experiments and of deterministic 
simulations (T = 0): (a) The roughness cr^(t) and (b) the 
variance A:^(t) of the local slope (only exp. 1 and simulations 
for clarity) as functions of time and fits based on the linearized 
TFEq. While cr^(t) can be fitted with the deterministic theory 
by adjusting the initial roughness ctq and the characteristic 
time scale to, this is not possible for the experimental k^(t). 
In the linear regime the deterministic /c^(t) is constant in time. 



mismatch in k'^(t) is most prominent for short times, i.e., 
for small cr^(t). In the following, we quantitatively anal- 
yse the early stages of dewetting in a linear approxima- 
tion of Eq. ([1]) to demonstrate the relevance of thermal 
noise in hydrodynamics on small length scales. 

In the beginning of the dewetting process the devia- 
tions ^(r, t) = h{T^ t) — ho from the initial film height ho 
are small as compared to ho- By expanding Eq. ([T]) in 
first order of dh and N we obtain a linear stochastic equa- 
tion in Fourier space. Its dispersion relation u;{q) = [1 — 
{q'^/Qo ~ l)^]/^o has its maximum at Qq = —^'\ho)/{2j) 
and a characteristic time scale to = 3/^/(7 /iq ^^q) (about 
300 s and 1700 s for exp. 1 and 2, respectively). The spec- 
trum reads C{q,t) = {5h{q,t) Sh\q,t)) = Co{q) 6^^^^)^ + 
^^W^ 4^ [e2^(^)^ - 1] with the initial spectrum Co{q) 
at t = 0. We assume the initial spectrum to be con- 



stant Co{q) = ^ctq foi" ^ < V^qo and zero other- 
wise. The spectrum has necessarily a microscopic cutoff 
qm = 27r/rm ^ qo at the scale r^, which is certainly 
larger than the fluid molecules. The cutoff for the mea- 
sured spectrum is set by the experimental resolution, e.g., 
by the pixel size or a smoothing length if the data is 
post-processed. Due to a rapid buildup of noise induced 
short wavelength roughness on a microscopic time scale 
tm = {Qo/Qm)^to^ the short wavelength part of the initial 
spectrum Co{q > qo) is irrelevant for the film evolution 
on the characteristic time scale to tm- 

In the case of a spinodally unstable film with ^''(/iq) < 
the dispersion relation uj{q) is negative for q > V2qo. 
This leads to two important differences between deter- 
ministic and stochastic dewetting. First, and largely 
independent of the initial roughness, for t ^ oo and 
q > V^qo one finds an exponentially decaying power 
spectrum C{q^t) Co{q) in the deterministic 

dynamics (T = 0) but we recover the algebraic capillary 
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FIG. 3: Slope variance from the experiments plotted ver- 
sus the roughness . The solid lines are fits based on the 
lineraized stochastic TFEq, whereas the dash-dotted lines are 
asymptotic fits according to Eq. (|2|. The influence of noise is 
reduced for smoothed as well as for less resolved data (exp. 2). 
The inset shows the same data in a double logarithmic plot 
highlighting the validity of the asymptotic Eq. (|2]). 
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for any finite temperature T. And second, for sufficiently 
smooth initial spectra which decay rapidly for q > V^qo^ 
the maximum of the deterministic spectrum stays at qo 
for all times, while the maximum of the stochastic spec- 
trum approaches qo from above as t ^ oo, see Ref. [5] for 
details. This noise generated coarsening process can last 
until non-linearities become important, effectively mask- 
ing the most typical feature of the linear deterministic 
regime, namely that the maximum of the power spec- 
trum stays at the fixed wave number qo. 

In order to further illustrate the spatio-temporal fea- 
tures of the dynamics we calculate the time evolution of 
the roughness of the film cr^(t) = / j^:^ C{q^t) shown 
in Fig. [SJa) for the experiments discussed above and the 
corresponding deterministic simulation. Both, the ex- 
perimental and the deterministic data can be fitted with 
the deterministic a^^Q{t) with parameters to and <Jo in 
reasonable agreement with the experimental parameters. 
The initial roughness ctq is experimentally not accessible. 
Thus, the time evolution of the roughness cr^(t) gives no 
clear indication for the importance of thermal fluctua- 
tions for the dewetting process. 

Yet, the variance of the local slope 27r a'^{t)k'^{t) = 
([V^(r,t)]^) = J^C{q,t) is sensitive to thermal 
noise. We find a strong dependence on the noise term 
for times t ^ to whereas the late stage k'^{t oo) = 
ko is given by the maximum of the dispersion relation 
ko = qo/{2 7r) alone. The reason why the noise term be- 
comes irrelevant for large t is the strong increase of cr^ (t) 
with time due to the deterministic growth of unstable 
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modes. Note that for a purely deterministic dynamics 
the variance k'^^^it) = changes in time only at micro- 
scopic time scales tm and remains constant in time if one 
chooses Qm = V^QO' In particular, the position q ^ k{t) 
of the maximum in the structure function C{q^ t) does not 
change during a deterministic dewetting process. In con- 
trast, thermal noise induces a time dependence of k'^{t). 
The experimental data cannot be fitted with the deter- 
ministic theory, because k'^{t) is clearly not constant in 
time. A simultaneous fit to the experimental data with 
ctq = yields kBT/{4:7Tj) and the cut-off wavevector 
q'^/qo as given in Table [D From the experimental pa- 
rameters we would expect ksT / {Ai:^) ^ 0.012 nm^ for 
both experiments, which is reasonably close to the fitted 
values. The values for q^^ correspond to cutoff lengths 
which are larger than the nominal pixel size. The rea- 
son is probably that the experimental cutoff is smooth 
rather than sharp. However, the ratios of the Vm for the 
different experiments matches quite well the ratios of the 
corresponding pixel sizes. In contrast to ctq, smoothing 
changes /c^ drastically. Deterministic simulation data is 
hardly affected by smoothing. Smoothing obviously re- 
moves fluctuations on small scales, which are exponen- 
tially fast damped by the deterministic dynamics but are 
permanently (re)generated by thermal noise in the ex- 
periments. A quantitative analysis of the experiments at 
different lateral resolution clearly confirms thermal noise 
as the most likely source for the observed /c^ dependence. 
The dependence of k'^ on the upper cut-off qm is a cru- 
cial test of the theoretical prediction as presented below 
inEq. 

As the microscopic cutoff qm is much larger than q^ 
there is a time regime t > tm ^ 10~^to up to t ^ to 
where k'^{t) is approximately given by 

k\t)^kl + J^^, with X=^|^, (2) 

independent of the initial conditions [5]. Since cr^(t) is 
monotonically increasing with time, the characteristic 
wavelength ^/27^/W(t) increases. Thus, thermal noise 
generates coarsening even in the linear regime for which 
the deterministic linear dynamics predicts a fixed char- 
acteristic wavevector. Fig. [3] shows k'^ as a function of 
cr^ fitted by Eq. ^ as well as by the full linear stochas- 
tic theory. Both fits give consistent values for k^ and x 
which are given in Table [D These values are also consis- 
tent with the values determined by the system parame- 
ters given above. 

From the fit values k^ and the theoretical relation 
q^ = —^"{ho)/{2'y) we get for the film thicknes ho in 
exp. 1 3.6 nm and 3.8 nm from the smoothed and raw 
data, respectively, and for exp. 2 /iq = 5.5 nm. These val- 
ues are in good agreement with our experimental data. 
Accordingly, one finds from the fit values for x the cut-off 
lengths Tm given in Table [J- in good agreement with the 
ratios of pixels sizes and smoothing length, respectively. 



It is remarkable, that one can determine reliable values 
for film thicknesses and lateral resolution by measuring 
the effect of thermal noise on fluid flow. This is a strong 
indication that the coarsening visible in k'^ {t) during the 
linear regime is indeed due to thermal fluctuations which 
can be described by the stochastic TFEq ([T]). Numerical 
solutions of Eq. ([Tj) would provide an excellent further 
test. 

In conclusion, the interplay of substrate potentials and 
thermal noise can strongly influence characteristic time 
scales of fluid flow and leads to a coarsening of typical 
lengths on a /am scale. The results will have impact on 
the future development of nanofluidics, enabling reliable 
predictions of fluid flow in confined geometries such as in 
liquid coatings or in lab-on-chip devices. 
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TABLE I: Experimental fit parameters for Figs. [2] and [3] 
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